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We use the variational approximation with double Gaussian type trial wave-functional approx- 
imation, in which we use the square root of the dispersion of the zero-mode wave-function as an 
order parameter, to study the out of equilibrium quantum dynamics of time-dependent second order 
phase transitions in (3-1-1) dimensions. We study the time evolution of symmetric states of scalar 
theory in several situations by properly treating the effect of the At/)* interaction. We also 
calculate the effective action and the effective potential of the theory with the precarious renor- 
malization. We show that the presence of a quenching of the mass-squared leads to second order 
phase transition nontrivially since the vacuum structure changes by absorbing the energy required 
for quenching, even though there is no symmetry breaking in the effective potential of the theory 
without quenching process. We also calculate the equal time correlation function, and then evaluate 
the correlation length as a function of the mass-squared. The time dependence of the correlation 
length varies depending on how the mass-squared changes in time. For constant mass-squared it 
gives the classical Cahn- Allen relation, and it leads to different relations for other time-dependence 
of the mass-squared. We also show that there exists a propagating spatial correlation after termi- 
nation of the phase transition process in addition to the correlation corresponding to the formation 
and growth of domains. 
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I. INTRODUCTION 

In recent years, there has been strong interest in the dynamics of the quantum fields out of equihbrium appearing 
in the evolution of the early universe or in the process of phase transition The formation of topological defects, 

the dynamics of inflation of the universe, and the inflationary reheating have been studied on the basis of scalar field 
model. The domain formation and growth Q were studied in the out of equilibrium second order phase transition 
using the Hartree-Fock approximation. During a second order phase transition, the time scale of relaxation of the 
scalar field lags behind the time scale of change of the effective potential. Consequently, the field evolves out of 
equilibrium as it tries to relax to a new vacuum, giving rise to a nonvanishing vacuum expectation value. Such 
nonequilibrium effects play a crucial role in topological defect formation in condensed matter systems as well as in the 
early universe ^ uSej^J ■ Kibble first showed how the correlation length is crucial in determining the initial density of 
topological defects Q ■ These ideas were elaborated by Zurek who proposed that it may be possible to quantitatively 
test the Kibble mechanism of defect formation in condensed matter systems such as superfluid He"* S|. He argued that, 
because of the phenomenon of critical slowing down near the phase transition point, the correlation length relevant 
for determining the initial density of the defects is not the equilibrium correlation length at the Ginzburg temperature 
but that at the time when the field dynamics essentially freezes out. He also found a power law behavior in the 
dependence of the correlation length on the quench rate. There have been various attempts to experimentally test 
the Kibble-Zurek prediction of initial defect density 9]. A comprehensive review on these issues is given in Ref. 1^. 

One crucial lesson provided by these studies is the importance of using time-dependent techniques to study processes 
of the system with explicit time-dependence. It has been repeatedly shown that classical and one-loop effective 
potentials are poorly defined and of little use in such dynamical systems. Another common theme is the importance 
of non-linear corrections to the linear dynamics. In these directions, the dynamics of spinodal decomposition in 
inflationary cosmology using the closed time path formalism of quantum field theory out of equilibrium combined with 
the non-perturbative Hartree approximation was analyzed in Ref. [TH . A feature common to all such phase transitions 
is the evolution far away from equilibrium and the exponential growth of the soft (long wavelength) modes, which 
necessarily lead to domain growth. Finite temperature field theory based on equilibrium or quasi-equilibrium methods 
does not describe all the processes of nonequilibrium evolution even though the imaginary part of the complex effective 
action yields the decay rate |l2|. To treat such nonequilibrium quantum evolution properly, Schwinger and Keldysh 
first introduced the closed-time formalism p^ . The closed-time formalism and the 1/iV expansion method have been 
applied to nonequilibrium A^"* theory to explain the phenomenon of domain growth [14i | . Recently, the Liouville-von 
Neumann approach has been developed that unifies both the functional Schrodinger equation for quantum evolution 
and the quantum Liouville equation for quantum statistical mechanics [T5l |. The renormalization of field theory in 
quenched second order phase transitions was discussed in Ref. 16] in this context. 

An obstacle in discussing the phase transition of the scalar field in (3-1-1) dimensions is that there is no symmetry 
breaking in the Gaussian effective potential of the theory with the precarious renormalization JJj, which is the 
only relevant renormalization scheme which does not require non-trivial assumptions such as the existence of a large 
momentum cutoff, large N limit of 0{N) scalar field, or the autonomous renormalization 0, 0|. On the other 
hand, there have been various studies on the 2nd order phase transitions of the scalar^ theory in (3-|-l)dimensions 
in relation to the OfiV) symmetric theory 0|, the non-equilibrium field theories [Tllll^. and the theory with high 
momentum cutoff [23. Most of these studies consider the region in which the phase transition begins, where the 
dynamical behavior of the unstable modes are the main interest. It is an interesting question to which vacuum the 
system will settle down after the transition since there is no vacuum with non-zero order parameter in the static theory 
described by the effective potential. In this paper we answer this question by showing that the vacuum structure of 
the scalar field theory changes when the second order phase transition occurs. 

The variational approach for the scalar field theory using the Gaussian effective potential was well studied in 
Refs. 10, I21I |23 |. and references therein. The renormalizability and the initial value problem for the Gaussian 
approximation of time-dependent scalar systems are checked in _23, 24] . Many authors have studied the symmetry 
breaking phase structures in the large N approximation psj ]. Non-equilibrium dynamics of symmetry breaking and 
the second order phase transition [13] have also been studied in the Gaussian frame work. The Hartree-Fock method 
has been popular and useful in studying the nonequilibrium evolution of the field jl4l ] . Even though the renormalization 
is well understood even for nonequilibrium fields in Refs. |14| a systematic and explicit renormalization scheme of 
the effective action has not been properly addressed for the systems undergoing phase transitions at least in (3-1-1) 
dimensions. In this paper, we propose such a renormalization scheme of the effective action for systems describing 
phase transitions. 

As field theory models for the second order phase transition, some exactly solvable models of the free scalar field 
were studied in Ref. |0], in which it was proposed that the exponential growth of the spinodal instability may ends 
at the spinodal line because of the back-reaction of the scalar field. The spinodal instability leads to the domain 
formation and growth which is determined by the equal time correlation functions. The formation and growth 
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of the domains were studied up to the point where the phase transition process ends. It is generally believed that 
the domain becomes static once it is formed without other dynamics involved. In this paper, by studying the spatial 
correlation functions determined by the approximation method proposed in Ref. ^ , we show that there exists another 
propagating correlations after the transition process is over in addition to the static domain. 

Recently, the present authors have developed a new quartic exponential type variational approximation (2S] which 
is suitable for double well type potentials in the quantum mechanical context. This approach was applied to the 
zero mode of the scalar field 27], where we have renormalized the equations of motion which describes the symmetry 
breaking phenomena starting from the symmetric states of the scalar theory in three and four dimensions. In the 
present paper, we use some of these results summarized in Sec. II. 

The article is organized as follows. In Sec. II, we briefly present the double Gaussian type wave-functional ap- 
proximation to the \(fA theory. Then, we introduce the Klein-Gordon like mode solutions, with which we rewrite 
the Hamiltonian of the system so that it includes the temperature of the initial equilibrium state. In Sec. Ill, we 
calculate the effective action and potential of the interacting scalar field with fixed mass in which the tachyonic mode 
described by negative mass-squared is included. It is shown that the typical effective potential develops a new local 
minima at the non- vanishing expectation value of the field at which point the mass vanishes even though the vacuum 
still presents at the symmetric point of order parameter. Sec. IV is devoted to the analysis of an exact mode solution 
of the field equation in which the mass-squared varies from positive to negative values. We obtain several asymptotic 
behaviors of the solution such as the instantaneous quenching, t oo limit, and ultra-violet (UV) limit. These are 
used in Sec. V to obtain the WKB approximation for the general field equation. In Sec. V, we study the self- 
interacting scalar field system with instantaneous quenching &tt — 0. After the quenching we allow the mass-squared 
freely evolve so that it increases due to the dynamical back reaction of interaction. We calculate the UV finite 
renormalized Hamiltonian and potential written in terms of a shifted order parameter, q, and mass-squared, m^, with 
an additional relation between q and m? . In Sec. VI, we introduce a large instability approximation, with which we 
analytically study the effective Hamiltonian and potential. It is shown that the vacuum structure of the theory is 
changed due to the presence of instantaneous quench at t = so that there occurs second order phase transition in 
the new effective potential. In Sec. VII, we calculate the evolution of the equal time correlation functions and present 
a simple formula which relates the correlation length and the tinic-depcndcnt mass-squared. We also show that there 
exists a propagating correlation in addition to the usual correlation which describes formation of domains. Summary 
of our results and some discussions are given in Sec. VIII, and three appendices are added at the end of the article. 



II. MODE SEPARATION OF THE SELF-INTERACTING SCALAR FIELD 

In this section, we summarize the double Gaussian type approximation of the self- interacting scalar field theory [23| , 
in which the back reaction of the field by the Ac/)'* interaction is considered in a natural fashion. Then, we briefly 
describe how the time evolution of the initial equilibrium state with inverse temperature (5 can be analyzed in the 
present framework. 

The Lagrangian of self interacting scalar theory in 3 -f- 1 dimensions is 



(1) 



where we explicitly included the volume integral so that we can write the volume factor in the zero mode part of 
the Lagrangian. In Ref. we have shown that the quantum mechanical excitation of the zero mode may play a 
non-trivial role in the symmetry breaking of the system, where we separated the zero mode from the other modes by 
the following form. 



J (fyi(t){x, t) 



-ijj{:>c,t) = 0(x,i) - (t>{t). 



(2) 



We use the unit ?i = 1 in this paper. Following the double Gaussian type approximation to the zero mode developed 
in Ref. jl^, we take the trial wave- functional 



^[(f>,ip] = A^exp 



1 



iTr{t) 



x(t) , , 



m(x,y;t) 



^(y) , (3) 
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where 



J (Px.. Then the Lagrangian leads to the effective action 
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n(k, <)G(k, t) - 2n2(k, i)G(k, t) 
^ ' G(k,i) 



(4) 



where D{y) is a smooth function of y defined in Ref. (22 



q{t) = y{t) - 



(</>2(t))2 



(5) 



describe the dispersion and shape of the whole double Gaussian type wave-functional of the zero mode, respectively, 

d^k 

.g»k (x-y)^j-]^^ t) denotes the width of the Gaussian trial wave- functional for the non-zero 



and (V'(x,i)V^(y,t)) ^ j ^^^^g. 
modes. The zero-mode potential is 



y>i, '7>o. 



(6) 



Note that the potential V is finite even in the limit of zero dispersion, q^ ^ or y 1, which is the distinct feature 
of the theories of infinite volume. The time-dependent variational equations are given by 



g+[m2(<)-7)2(i)]g + 
G{k,t) = ^G-\k,t) 



^{y-3)q^ = 0, 

^^G-\k,t)G'{k,t) 



where 



2[k^ + m^{t)]G{k,t), 



G(k,t). 



(7) 
(8) 
(9) 

(10) 



A distinctive feature of this equation is that the square root of the dispersion of the zero mode, q, satisfies the deformed 
classical equation of motion ©, and can play the role of order parameter. 

The potential divergences in the integral G(k, t) of the equations of motion are absorbed into the bare mass 
and the bare coupling A leading to the renormalized ones. We have shown in the previous paper QJ] that 
the equations of motion are renormalizable. But the proof is not complete in the following sense: 1) the theory 
considered in the previous paper does not include the unstable modes where negative mass-squared appears and 
spinodal instability grows. 2) the proof does not include the renormalizability of the effective action even though 
it shows the renormalizability of equations of motion. In computing the effective action, we have extra divergent 
integrals of the form J^^G^^ and Jj^ aj^G(k, t) which may cause extra complications. For simplicity we set y = 1 in 
this paper so that it is non-dynamical. The dynamics of y may give non-trivial contributions to the effective action, 
however, its inclusion is not difficult because dynamical equation (jT)) for y simply relates y with l/<?^ in the precarious 
renormalization scheme. 

To consider the spinodal instability, we rewrite the equation for Gk by introducing the mode solution (/'k(i) of the 
Klein-Gordon equation. 



(11) 



where (f>i^ = i?ke*^'' With the identifications, 



^k^k 



1 

2' 

k^+^l^{t) + ^q'{t) + 



A 



G, 



(12) 



we obtain 



G(k,t) = ^l{t)Mt)- 



(13) 
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For time-independent system with fixed mass m(t) = > 0, the solution for 0k(^) is given by 

1 



•\/2Wi,k 



(14) 



where Wi^k — y'^^ + mf is the initial frequency of mode 0k- The direct link of the mode solution and the Gaussian 
approximation was discussed in Ref. [2^ and this Klein-Gordon type mode solution was discussed in RefJ3 in the 
context of Liouvillc-von Neumann approach to the scalar field theory. As in the case of the reference we can 
define the thermal expectation value at t = 0, which gives the Hamiltonian density 



H{q,p;G,G) 



1 



A 



2 ' 2"^ ' ' 4!'' 
where f3 is the inverse temperature at t = and 



i/k(i) coth 



A 



<Pk{'t)<P^i't) coth 



f3uji. 



A 



Luiit) = y+ f/it) + -q\t) + - / 0Ut)0k(i) coth 



(15) 

(16) 
(17) 



The generalization to thermal state comes from the equivalence of the Liouvillc-von Neumann equation and the 
Gaussian wave- functional approach, which was well described in Ref . 1151 . We do not digress about this point in this 
article. This Hamiltonian can be compared with Eq. (7.18) in Ref. [I^- Note that the frequency in the presence of 
temperature includes the temperature dependence in the momentum integral. 



III. EFFECTIVE ACTION FOR THE CASE WITH CONSTANT MASS-SQUARED INCLUDING THE 

TACHYONIC MODES 

Since solving Eqs. JHJ, and (|11() is not simple, we start from the simplest case, the one with the mass-squared 
m^(t) = m? being fixed irrespective of its sign. Even though, in realistic system, the mass-squared may change by the 
interactions, we assume that the mass-squared is kept fixed in this section by an external mechanism which we 
do not specify. This procedure may need an external work, which may alter the energy of the system. For simplicity 
of the calculation, we calculate only the zero temperature case. In the subsequent sections, we calculate the case of 
time- varying mass-squared ignoring the self interaction, and then we synthesize the two cases to calculate the effective 
action for a self-interacting scalar field theory with time-dependent symmetry breaking. 

We first define an integral notation, /k(i), by 



/ /k(t) 



We always use the notation: rn? — \tv?\ if m? is negative. During the calculation of the effective action we frequently 
encounter integrals of the form: 

lNim')^U {k'+mT-'^'. (19) 

where we keep the square in the argument since the mass-squared can be negative. In the case of positive definite mass- 
squared, > 0, this integrals are well studied by Stevenson 17] in which the reduction formula of the divergences 
are given by 

/o(to^) - Io{mji) 



|k|=m (27r)3 



hit), -m^{t) = fh^it)>0, 



=0(2^ 



(18) 



m'^{t) > 0. 



= 2("^' - ^IW^l) - ^i^' - mlfl-iiml) + ^i3(x), (20) 
= -^(™'-"^fl)^-i("^fl) + ^i2(x), 
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where 



and 



L2{x) ^x\nx - x + 1, Lsix) = -[2a;^lnj: - 2{x - 1) - 3(a; - 1)^]. 



(21) 



Beside the above formula we need the foUowing two formula which relates the integrals with positive and negative 
mass-squared, 



/o(m2) + im2/_i(m2) 



167r2' 



(22) 



Armed with these formula H2U|I and H22|l . we are prepared to calculate the effective action. The Klein-Gordon 
equation Hll() can easily be solved to give 



where uj^ 



for < -m^ > 



1 



1 

/2wk 



[cike^-* + cake^^-*] , < k^ < -m^ 



(23) 



-^u^^t^ k^ > -m^ > 0, 



k^ -I- m^, and cik and C2,k are arbitrary constants of 0(|k|°). For the 



stable modes, k^ > — m^, only the positive frequency modes are chosen to give time-independent Gk, 



Gk 



[Icikl'e^^--* + Icakl'e-^^-* + cikC^k + c*ikC2k] , < k^ < 

2wk 



k^ > -m^ 



(24) 



In this paper, we follow the precarious renormalization scheme |l7l ITgl |2^ . which uses the coupling constant 
renormalization condition, 



1 1. / 2n 1 



From this, the definition of mass in Eq. H10|l becomes 



where we explicitly denote the t dependence in m^{t) for later use and 

f = 9' + ^Gk-/o(m2) 



(25) 



(26) 



(27) 



is a finite redefinition of variable q, which plays a crucial role in the renormalization of the effective action. In the 
present simple example, = q^ + ^{~^) /|k|<m^k, where q = q for positive mass-squared. Later in this paper, we 
write the effective action and potential in terms of q since it is simpler. Subtracting the equation at the renormalization 
point, {qR,mfj > 0), from Eq. (|26() and using the integral reduction formula (|20|l and (|22l) we get 



'z'-9fl + ^(a^lnk|-a^ + l) 



(28) 



This equation can be simplified by introducing dimensionless variables, <f> = 167r 



(1 - k){x-1) - ($2 _ $2^) ^ xIti\x\. 



2 Q 



and K = — 



327r2 



(29) 



The equation l|29|l can be solved by graphical method, i.e., by explicitly plotting each side of the equation as a function 
of X. As an explicit example, we consider n > I case in Fig. 1. The equation has three roots for —k + 1 + e~'^ < 



7 




FIG. 1: Graphical solution for Eq. The thick curve represents the function a; In [a;|, the thin curve is the effective potential 

V at a given time, and the straight lines A, B, C represent the left hand side of 1291 1 with — <l>fj = 0, k — 1 — e~'^, ft — 1 + e~" 
respectively, for k > 1. The arrow indicates how the lines move as increases. The potential has a local minimum at 

a: = 0. 



<i>2-<i>|j < -{-K+l + e^'^), two roots for <i>^-<i>|j = ±(-k+ 1 + 6^**), and one root otherwise. For 
the straight hne passes through the origin. Note also that <i>^ — $^ should be greater than —^j^- The improvements 
in this solution from the previous works are the following: First, we always have roots of Eq. H29|l for any value of 
> 0. Second, a: = is not a special point now since the potential is continued for negative x. 

The divergent part of the effective Hamiltonian H15() consists of two parts: l^^-^l^ 's i^J '^^'^^ -^k ^gm- 

is convenient to change the first part using Eq. 1)27(1 to give 

W -U [gS^ + ^(<?^ - - ^lUm') + - f), (30) 



2 8 VA y 2 " 2 " ' 8 ' 2 

where the divergences in /q and bare mass, /i^, are separated from the finite contributions. If one incorporates the 
renormalization condition (|25|) . the last term in Eq. (|30|l goes to zero since the bare coupling is —0. In the second 
part, we separate the time-dependent part from time-independent part using the explicit formula for Gk, Eq. 124|) . 



k 2 k 



where the time-dependent part is 

"1 

'|k|<m J|k|<m 

Ignoring terms which go to zero when A ^ — 0, we get the effective Hamiltonian 



J\k\<fh J\k 



+ (cikC2k + C^kC2k 



Gk'- (32) 



Heff{q,p) - ^ + ^{q'-f) + ^f+h{rn')~^lUni') + e{-x)m%x'VL{fh,t) (33) 
= y + ^(g' - f) + 0(-x)m%x^VL{rn,t) + VD{q,m'), 
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where the divergent part, Vd, becomes 



(34) 



D + ^[f ~ fn] + him') - h{m%) - ^[/qK) - Io{ml)] 



327r2 

.4 



D + ^mUf - ql) + ^L,ix) - ^(x - ir 



D + \mlx(f - ql) + ^[L3{-x) + 2x]-^ix- 1) 



2 

X = — 5" > , 

x= < 0, 



with D — + — -3-^0 (^fl) being a 5 independent divergent constant. Let us define the normahzed effective 

2 8 
Hamiltonian by 

Heff - D 



n 



m|j/(647r2) 

T^ + Oi-x) 

327rV 



-2^, .^.2x._ .^ , -i^JLP + ££^3. _ 2x<i>2^ + 2L3{\x\) + 4x0{-x) - k{x - l)^ (35) 



= 6ATT'9{-x)x'^VL{fh,t) + 



+ ei-x) 



-R 



Mii'x^VL{m,t) 



MTT^x^VL{m,t) 



■yiTx'x 



R J\V\<m 



■yiTx'x 



m 



R J\\A< 



+ x($2-ci>|)-i(a,-l)(x-l + 2/t) 



x2ln|x|+(^-K)(a;2-l), 



where we write the Hamiltonian in several different forms using Eq. (|29|l . 

Let us consider the case where the initial value is given by C2k — and cik = ci , to present an explicit form of the 
effective potential. Then, VL(m,t) exponentially decreases to zero as time t — > 00 and 



|k|<»7i 



|ci|2 |ci|2 
Gk = — / r-^^ dk = — 



47r2 7q ^m'^ - k? 



(36) 



where the function / is maximized at = — , and the integral can be performed using the steepest descent method 
to give (See Appendix A): 



|Ci| /TTTO 2mt-l/2 



(37) 



|k|<m 



(38) 



47r2 V 2^3 

Therefore, the renormalized effective potential at a given large time t becomes 

V(a:) = 7^|p^o = 8|ci|V5(mfl(-x)i/2O0(-x)-a:2ln|x|+ Q-K^ (x^-l), 

where g{y) — 4 /— ^^e^^^^^^. The renormalized potential is symmetric about a; = if ci = 0. The first term of the 
V V 

effective potential (|38ll exponentially increases as time so that one cannot have a large negative mass-squared at a 
later time. 

Now we discuss the shape of the effective potential as a function of q as seen in Fig. 2. Since g = g for m2 > 0, 
this gives the effective potential for q for positive mass-squared case. As can be seen in Fig. 1, there is at least one 
root X for a given q. In the case of two or three a: roots, the true x values should be chosen such that it minimize 
the effective potential (|38|) . Then the absolute minimum of the potential is at $2 = q. The potential increases as a 
function of $2 until (f>2 — $|, — e^" — 1 + k where x = e~'^ . It then decreases until <I>2 — <I>^ = k — 1 where a; = 0, 
which is a local minima of the effective potential, and then it increases again. The metastable state corresponds to 
a; = so the theory is massless at that point. The effect of the inclusion of the negative mass-squared modes becomes 
apparent for large q' , in which, the potential exponentially increases in time and the mass-squared is negative. This 
means, for a realistic system, that the modes with negative mass squared soon decays to positive stable modes since 
it needs infinite energy to keep the negative mode continually. 

Some comments are in order. Due to the inclusion of negative modes the solutions to the equation H29|) always 
exist. Stevenson [l7j | found that the present technique is inadequate to find the effective potential for negative k for 



V 
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FIG. 2: Typical form of the effective potential V at a given time t. The horizontal axis represents q. The global minimum at 
the origin is at g = and the local minima at g 7^ is the point of q which corresponds to x — 0. 



the following reason: At x = 1, we should have (m^,g^) = {rn\,q\). There are two solution for x to the Eq. H29|) if 
K < and the solution at a; = 1 corresponds to the larger potential value of the two. This does not obey the guiding 
principle which chooses the x which gives lower potential value. To overcome this difficulty he used a rescaling of 
variables which maps k < into < k < 1. This is also the case for the example above. 



IV. SPINODAL INSTABILITY AND FINITE SMOOTH QUENCH BY CONTROLLING THE 

MASS-SQUARED 

In this section, we consider a finite smooth quench model, which was previously considered in Ref. |15| . The finite 
quench transition is described by a scalar field with the mass given by 

m| - m f mf + TO f f t\ , ^ 

m\t) = ^ ^ - ^ ^ tanh j , (39) 

where to| and TOj are both positive definite. At earlier time, t = —00, the mass-squared has the initial value m| and 
at later time, t = 00, the final value —iTif- Here r measures the quench rate: the large r implies that the mass changes 
slowly, whereas the small r implies a rapid change of the mass-squared. In particular, the t —^ limit corresponds 
to the instantaneous change from m| to —mj at t = 0. To find the Fock space for each mode one needs to solve the 
classical equation of motion 



„ mf — TO f mf + TO f / t 

c + — —: — - tanh 



Mt) = 0. (40) 



It should be noted that the long wavelength modes, |k| < my, lead to the change in the sign of the frequency at later 
times {t ^ t), 

ulit) = k2 - m} < 0, (41) 

and suffer from spinodal instability. Each long wavelength mode has a different quench time determined by w^(tk) = 0. 
The solution to Eq. (|40() are found separately for stable modes and unstable modes. The stable modes {k > mf) have 
the solution 

0k(O = — ^e-'"--* 2-Fi(-«™+,k, -«™-.k; 1 - «™».k; -e'*/"), (42) 



where 



W±,k = TT^' (43) 
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with 

uj,M - \/k2 + m2, t^^-k = - m2. (44) 
Whereas the unstable modes (|k| < m/) have the solution 

<^kW = — ^e-^'^->* 2i^i(™k, -T<; 1 - z™,.k;-e'*/"), (45) 



where the complex frequency is given by 



i«k = ^^^-^ -■ (46) 



At earlier time [t ^ — r) before the quench begins, both solutions (|42|l and H45|l have the same asymptotic form 

0.,k(t) = -^e— (47) 

V ^^i,k 

which matches the initial condition H14|l . At later time {t 3> r) after the completion of quench, we find, by using the 
asymptotic form of the hypergeometric function 30], the asymptotic form for the stable modes 1)42(1 

hsMit) - C^S+,k^== + C5-.k^==, (48) 



where 



J^TO-u^^^^^ (49) 



and the asymptotic form for the unstable modes H45|l 



huMit) - Ca+,k^== + Cc;_,k^==, (50) 



where 



^ /^/,k r(i - zLj,,kT)r(^/,kT) /t^/,k r(i ~ ^u^^.kT)^(-tJ/,k^-) 

y t^i,k ^«k'^r^(WkT) y w^^k -wkTr-^^-u^k-r) 

From these results, the transition rate from the initial vacuum to the final vacuum and the coefficient of the 
negative frequency solution leading to the particle production rate were calculated in Ref. [l^ . The asymptotic form 
for 4'*fu k'^fUM was also calculated to obtain the scaling relation of the domain size 




^D{t)^\l—. (52) 
where 

< = i-y[C(3)-l]|«;kp. (53) 

For use in the next section, we explicitly present some more results not mentioned in Ref. The asymptotic 

form of the absolute square of the stable modes are 



G/.,k = 0/„k(i)</'/5,k(i) = 2^ [1 + K cos (2a.;,ki - ^^)] , (54) 
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with the time-independent constants g^, R^, and tp^ given by 



si 



sinh^ TTW+.kT + sinh^ TTU^^kT 
sinh TTOJi.kT sinh ttlo f_]^T 



sinh 



sinh 



In 



sinhTTW-.kT sinh TTOj+.kT 
r(l — ia;/^k''')r(l + ja;+^kT)r(l — ia;__kT) 



r(l + ia;/^k''')r(l — iw+,kT)r(l + iw_,kT) ' 



(55) 
(56) 
(57) 



where i/'k ^ for |k| ^ l/r, and — 1, = for |k| ^ l/r. On the other hand, the uhra violet hmit of 
'^/s,k(0'^/s,k(0 is given by 



2uj 



1 0(6-""-"^), k2 > ^ 



(58) 



Using the properties of F functions we get the instantaneous quench hmit (r 0) of the stable modes and the 
unstable modes: 



lim 0/s,k(O 

T— >0 



lim k(0 

T— >0 



1 



1 



t^i.k 



g-iw/,kt _|_ 



^/,k 
^^i.k 



^/,k 



and their absolute squares are 



l^i.k 
Wi,k 



(59) 
(60) 



limG^s,k(i) 
limG/,,^k(i) = 



where 



±,k — o 



1 I t^/.k _|_ ^i,k 



2 V w. 



(fi+ k — f^-.k COS 2a;/_ki) , 
f2_^k + ^+,k cosh2aj/^ki 



1 f W/,k _|_ ^i,k 



2 V w. 



(61) 
(62) 

(63) 



All of these limiting properties are used in the next section to obtain the WKB mode solutions with arbitrary time- 
dependent mass. 



V. TIME EVOLUTION AND RENORMALIZATION OF A SYSTEM AFTER INSTANTANEOUS 

QUENCHING 

The model described in Sec IV does not describe a real system undergoing second order phase transition because 
the spinodal instability increases indefinitely. This model describes an intermediate process of the realistic phase 
transition toward the spinodal line, m'^{tk) = k^, at which point the instability of the mode (/)k stops to increase. In 
this section we consider a self-interacting scalar field model by properly treating the back-reaction of the field through 
the X(p^ interaction. All the instabilities end to increase at time T defined by ra^iT) = 0. Let the mass-squared of the 
scalar field m?{t) be positive value mf for t <Q and then it suddenly jumps to negative value —rn^f at t = t <^ ^i^- 
We consider the quench to occur in time r which is far smaller than the normal time scale (1/m), which is why we 
use the word "instantaneous" even though it is not zero. And the time r should not be zero since no realistic physical 
processes can alter the mass-squared instantaneously. Unlike the example in the previous section, we do not keep the 
mass-squared, m^{t), fixed, instead we let the mass-squared vary freely by self interaction after the quench t > t. The 
self-interaction of the scalar field leads to the gradually increasing mass-squared to a positive value as time. 

The equation of motion for the mode solution (|11|) . then, becomes 



(64) 
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Since the mass-squared is negative at t = r, the modes for t > t are divided into two categories depending on the size 
of momentum relative to the mass 



c(0 



4, k2 > -m^{t). 



(65) 



In the WKB approximation, the unstable mode, 0^, has exponential type solution and the stable mode, <f>^, has 
oscillatory solution. Let us assume that = + m^(t) increases for t > t due to the back reactions. Then, 

the unstable mode, 0^, becomes stable at time ife, given by k^ + m?{tk) — 0. Conversely, for a given time t the 
modes 0^ with |kp < —m'^{t) are unstable and the modes (j)^ with |kp > —■m?{t) are stable. Since the mass-squared 
increases as time, the number of unstable modes decreases and finally the unstable modes disappear at time T, where 
m?(T) = 0, which means the ending of the instability growth. Now, the dynamical process can be divided into three 
cases depending on the regions of time. First is the time before the transition i < which is given by the initial 
condition, second is the spinodal development process, < t < T, the phase transition, and the final one is the 
stabilization process after the termination of the instability growth, t > T . 

Let us assume q{t) = and w?{t) = > for < < 0. Then, the initial state becomes 



and Eq. ifT^ yields the gap equation 



Mi < 0) 



2uj, 



1 f P^iM 

■ coth 



(66) 



(67) 



where is the initial bare mass and A is the bare coupling. The infinities in the k integral are absorbed into the bare 
mass and the bare coupling leading to finite and positive mf. For convenience we use the notation Gj^it) = — cij^(i) > 
for the modes k^ < — m^(t); 



V(t)_k2-^^0*(i)0k(O coth 



> 0. 



Using Eq. H59|l and the WKB approximation of Eq. H64|) . we get the solution for a given time t > t, 



lit) 



where we use the matching condition. 



Cs^. 



1 



2^Ai[-^fe(t- ife)] 



< k^ < -m^{t), 
-7ii?{t) < k^ < 
\s? > ml. 



(68) 

(69) 
(70) 



1 f'fc 
— pjt 



^{t')dt' 



1/2 

TT Bi[-/Ufc(t - tk)] 



1/2 



+ 

(7,k 



Wk(t) 

1 



V2^k(i) 



: sm 



Mt')dt' - - 



LUk{t')dt' - - 



(71) 



of WKB approximation at time t = t^., where /i^ = 
in this paper, we use the notation 



(tk) 



1/3 



, and Ai(a;) and Bi(a;) are the Airy functions. Later 



o^{t)= / Mt')dt', e+{t)= / Mt')dt', 9^{t)= / tokit')dt'. 

Jo Jo Jtk 

In the case of the unstable modes, |k| < m(i), the frequency satisfies |k|r <C 1, which leads to 

-ek(t) 



it) 



1 



^i,k 



l^i.k 
'^/,k 



W,;.k 



(72) 



(73) 
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Therefore the WKB extended stable mode for —m?{t) < < rni becomes 



1 



I^i.k V ^/.k 

Similarly, the stable modes with moderate frequencies (my < |k| ^ 1/t) become 

1 



0k w 



2V2wk(i) 



^/,k 
t^j.k 



^»,k 

^/,k 



^/,k 

t^i,k 



^/\k 



and the stable modes with very high frequencies (|k| ^ 1/t) become 



^i.k 



^+.k 



V2wk(i) 



(74) 



(75) 



(76) 



As shown in Eqs. (|75|l and (|76|) the UV behavior of the stable modes takes different form depending on the value of 
the large momentum compared to l/r. 

The absolute square of this mode solution, Gk(i) ~ <Pt:{t)'P'k{t), can be written as 



Cl^^k + f]+,kCosh26'k(i) 



,k^< 



.9k 



2wk(i) 

9i 



where g^, R^, and 9^ are given in Eq. H55|) . and 



[l + cos (261^ (t) - )] , -m2(<) <k^ <m} 
[l + i?+cos(20+(t)-V^)] , k2>m2. 



ffk 



,2ek(tfc) 



'2ek(tfc) 



i?,7 = 



tan — 



16 



,k 



4e2«k(tfe) -I- e-2ek(tfc) 



1/2 



.k 



4 



Note also that G^{t) in the UV hmit (|k| > 1/t) is 



1 



2a;k(i) 

On the other hand, in the intermediate region (nif < |k| ^ l/r), it becomes 

1 



Note also that the time derivatives of G's satisfy 



[r2+,k + f^-,kCOs26l,[(t)] , m/ < |k| < l/r. 



(77) 
(78) 



(79) 
(80) 
(81) 

(82) 
(83) 



G^W = TT^G^ (i)+f^+.ksinh20k(t), 
^^k 



Therefore, the i/k in Eq- 116(1 becomes 



m2(t) 



2-.G^(t)-g±i?± sin [20±(i)-V'^]- 



2c.2(t) 



32^4 



+ £7.9k^i?J Bin [2e^it) ^J] + i.9±^k + ^ [1 - 5^(1 - )] (Gf) 



32^^"^ + ^f^+,ksinh20k(<) +ci^Gl^ - -(ikl^-,k. 



(84) 
(85) 

(86) 
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where we can ignore the last term of H86|l for > rnj since (1 — 



sinh^ TTW-)- ^ sinh^ 7rti;_,k'''l 



sinh TTcj ^.kT sinh 7rcj,;.kT 



1 for 



most range of |k| > m. 

The ultra-violet divergences are related only to (j)^ modes. From the structure of written in Eqs. 1)82(1 we see 
that the only divergence in Gj^ comes from l/[2a;k(0] term. It is already shown that the equation of motion for this 
form of Gk is UV renormahzable 0) 113 ^^^^ the condition given in (^5)1 . Using Eqs. (|Sn)l and we get the 

effective Hamiltonian 



H^Hj + + him') - ^I^im) + \ 



l)wk 



(87) 



where 



Gk 

32cj4 



/ ^sin[2^+(0-^k"]+/ 



^^sinh24(t) 



k2<- 



^k<-^k 



(88) 

'^k^^-,k 



k2<- 



The ± sign in , R^^ , and 9-^ should be chosen by the sign of |k| — mf. The large momentum expansion oi — 1 
becomes 



' (a;/,k - ^i,k)^ 
2a;i^k^^/,k 
g+-l= i (mf+m2)2 



8c^4(i) 



m < |k| < 1/r, 

m < |k| < l/r, 
|k| » l/r , 



(89) 



which gives a large contribution proportional to ln(m/jr) to the integral. This logarithmic contribution to the energy, 
which is related to the "instantaneous" quench process, changes the vacuum structure of the system in the sense that 
the transition probability from the initial vacuum to the final vacuum goes to zero in the r ^ limit. Since such 
zero-probability transition is not physical, we should restrict r 7^ so that the transition probability between the two 
vacuum does not vanish. Let us define the integral. 



/_i(m^,TO/) 



,1/t 




Jnif 




1 

4^ 


— In m f 


1 

4^ 


— In m f 



+ \Al^f, 



m \ 

1 + ^ + In 

mf 



1 + Vl + m'^T 



1 + y^l + m^ /m 



1/2 _ l + ^l + \x\/xf 



(90) 



where Xf = m'j/mj^, we ignored terms which vanish in the t limit in the last equality, and we defined the 
integration range from mf to avoid the infra-red (IR) divergence which appear at the k with u>-k{t) = if m'{t) < 0. 
We denote the finite part of the last integral in (|87|l by 



Vr 



ig+ - l)wk(i) 



+ ,k 



(mf+m2)2 

-^/_i(m^m/). 



(91) 



Then, the Hamiltonian becomes 



l^[m) - -I^{m) + ^^I_^{m\mf), 



(92) 
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where Hj + is the UV finite part of the Hamiltonian. Since the form of the divergences of the Hamiltonian is the 
same as Eq. (|33f) . the rcnormaUzation can similarly be done: 



H = D + Hf + Vr 



647r^ 



a;(<i>2-ci>2j)--(a;-l)(a;-l + 2«) 



(mf 



(93) 



where the divergent constant, D = -;:;IJ'^qR + Iiim^) + -^/q (rn/?). The last term in this Hamiltonian contains a large 

2 8 ^ 

constant contribution, which is the energy added by the instantaneous quench at t = 0. 

If an observer is confined for t > t and the experiments is done with the energy scale ^ the experiments 

cannot probe the existence of the very high frequency behavior (|82(l of Gk- In this case, one may do additional 
renormalization of the InTO^r term which is related to the sudden change of the mass-squared at t = 0, 



H 



D' + Hf + Vr 



6Att' 



where 



1 + 2k) + 2{xi +xff 



-l + (l + N/x/)-i/^-ln 



1+^1+\X\/Xf 



(94) 



(95) 



j /m 



= m'l/mjj, and D' = D — 



(mf 



4)' 



■In 



r. The divergence in /_i(m^,TO^) is the origin 
In some sense, Vvac denotes the vacuum 



327r2 

of the time-dependent renormalization, Eqs. (5.8) and (5.9) of Rcf. 
structure since it comes from the UV behaviors of the mode solutions which will not be affected by small excitations 
of moderate modes. In this point of view, Hf + V^- may be considered as an excitation to the vacuum state. It may 
be helpful to figure out the structure of the potential V^ac- The new term Vyac — V always decreases in x, which leads 
to the possibility of symmetry breaking. Since the mass-squared satisfies Eq. H26|) , the equation which determines x 
by $ is given by Eq. (|29|l and its solution is obtained by graphical method in Fig. 1. 
Note also that in Eq. (|27|l the difference, 



0i-x) 



|k|<m 



1 

2tJk 



[5^-l + .9j^^cos(20±(i)-V^)] 



(96) 



is positive definite on time average due to Eq. (|89() . Therefore, the minimum value of should be a positive number. 
Note that the integral /^Gk/i^k 1^ is apparently IR divergent at the value of k where k^ -I- = if 

< 0. As one may see in the previous literatures [23. |27| . there is no IR divergence in the equations of motion 
expressed in q and Gk. Therefore, this IR divergences come from the approximation of Gk using the WKB solution. 
This origin of the IR divergence reminds us the validity range of the WKB approximation. 



fvw?{t) 



«1, 



(97) 



which clearly fails to be satisfied at "the classical point" t = tk where Wk(^fc) ^ 0. In this sense, the IR divergence 
is the artifact of the WKB approximation around the transition point of ij^klO- The IR divergence does not have 
physical origin but comes from the bad choice of the solution of (f>k near t = tk- As an evidence of this argument, if 
one uses the exact solutions (Airy functions) near t = tk-, no divergences appear. To remedy this IR divergence, we 
need a generalized WKB approximation |^, which is beyond the subject of this paper. 

Fortunately and in fact naturally because this IR divergence comes from the failure of the WKB solution to 
satisfy the validity range of the WKB approximation (|97|l . this IR divergences present only in the kinetic terms for 
m^{t) < 0, which do not influence the effective potential at a given time. The equations of motion which determine 
the Gk functions are also IR finite, which is related the physical correlation length. Moreover, the effective action 
after the termination of the spinodal transition, t > T, is also IR finite. In the present paper we restrict our interest 
only to these cases. 



VI. LARGE INSTABILITY APPROXIMATION 



In the last section, we renormalized the effective action and potential. In this section, we calculate the renormalized 
effective Hamiltonian and potential for the case T 3> 1/m. The validity of this approximation should be checked by 
the dynamics at the time t < t < T, which needs the full understanding of the IR properties mentioned at the end of 
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the last section. In the absence of this knowledge, we can simply assume m/ to be large and the renormalized coupling 
to be small so that it takes long enough time to increase the mass-squared to zero from —rnj:- This approximation 
leads to the unstable modes to have large instability because of the exponential increase of the WKB solution, which 
is the reason we call it the large instability approximation. 

Most part of this section is devoted to the evolution of the system at time t > T except for the short calculation of 
the effective potential for < at the end of the this section. Therefore, the mass-squared is always positive definite 
and there is no unstable modes. The function, in this approximation, is dominated by the exponentially 

growing term given by 



and its integral over k. 



|k|<m2 



Gl{t) 



(98) 



(99) 



is approximately integrated by using the steepest descent method in appendix B. The resulting integral formula are 
summarized in Eqs. HB7|I and HB11|I . Before proceeding further, we define some parameters 



\2TiatJ 

- (^)". 

K = k + ikQ{t), 



(100) 

(101) 
(102) 



where I is a large scale which determines the instability, 1 < a < 2 is a number, and t is a time scale determined by 
fc^ -I- 'm?{t) — 0, both of which arc implicitly dependent on k and Q{t) = Q(t, r = 0) defined in Eq. IjBlOjl . 
We can approximate the integral in Eq. H99|l as 



+ ,k 



|k|<my 



-2e^(t) 



2u;«(0 



2w„.(0 



(103) 



Therefore, 



,fe 



-2e^(t) 



n 



+^g2jer(t) 



-t- Q{m,t), 



(104) 



,2^fe(i) 4i 

where Q{m, t) ~ G^{t) — Io{m'^), of which we present a rough estimation in appendix C. We do not need the 

full expression of Q{m,t), rather it is enough to know that it behaves as Q{m,t) cx — mj-y/x/x^ for large x. 
Other integrals in Eq. (|88|l are given by 



cm 



I ■ 



-2e^(t) 



2i 



^ + ,K 



n 



2i^l{t) 2^l,{t) 



2u;m ^2col,it) 



where in the second integral we use 

V'k ^ 7r/2, 5k - il+,ke2^''(*'=), and i?^ ^ 1- We ignore the integrals 



and 



' |k|>m 



ffk -^k 



|k|>, 



(105) 

(106) 
(107) 

cm 



sin [29^^ (t) — tp^^ ] compared to its small momentum term in Eq. p05f) . 
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Therefore, the Hf becomes 



Hf 




(108) 



'2»e, (t) 



-2i9r{t) 



From Eq. (|89|l we know that the last integral is finite and small compared to the exponentiated terms. Similarly, the 



integral /, 



|k|>mj 



can also be ignored compared to other terms. 



We approximate Vr also by its value around k — irif. We should also calculate - / 3k in Eq- H87|l . which gives 

2 J+.k 

1, (109) 



X 











--) 











X 



ijg'{xi,Xf) HE 
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where g' is a positive constant and the integrals needed in this equation is calculated in Appendix C. 



One can sum everything to obtain i?e/ / — D' + Hf + Vr' 



647r2 



Vvac- Before doing this summation, let us investigate 



where does the main contributions come from. Since loj. ~ rn{t), it is enough to observe the m(t) dependence of each 
terms. The main contribution comes from the unstable modes which have I factor if the value of x is not very large. 
For larger value of x the —x'^ In \x\ term in V^ac starts to compete with I term and then dominates the potential for 
larger x. Therefore, to know the behavior of the effective Hamiltonian, we need to keep these terms only: 



eff 



D' 



K 



m 



+ .k 



2^ 



(110) 



8i 



n 



±:^g2»er(t) 



2w„(t) 



647r2 



rX^ Inx, 



where K represents the sum of all kinetic terms and oj^/mp, — \/x + x with x = k"^ jm^. Therefore, the effective 
potential for very large t and a; > becomes 



V = 



eff 



D' 



my(647r2) 



2v , , 

— \\/x + x 



2VF 



— X Inx, 



(111) 



where v 



2^-KH[Xi + X/) 



and x ^ 1 for large t. Note that a natural stabilization of the effective potential occurs 



for large t > T due to the factor e"^^ (*) 

Until now, we did not calculate the effective action for negative x which is impossible until we treat the IR divergences 
properly, but the effective potential can still be calculable which exponentially increases in time as can be seen in Sec. 
III. The main contribution of the unstable modes comes from the integral Jj^ 



ik2<-m2 <^k^k which becomes 



CJl(;<-^k 



— / dfctik(t)f)+.kfc'e2«~''(*) 



(112) 



,kt 



TTk? 

27r 



3/2 



where kt is determined by kj = / Due to the exponential factor in Eq. (|112|l . the potential 111211 for 

large t is very sharply inclined to the vertical axis for negative x so that x increases. In summary, for a given time i, 
the effective potential for x has minimum at x = and increases in both directions. For large x ~ u^/-^^ the x^XvlX 
term starts to compete with v^fx. For larger x > u^/'^, the potential starts to decrease. The relation between _ 
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FIG. 3: Graphical solution for Eq. The thick curve represents the function a;ln|x| and the thin curve represents the 

effective potential V/v (we set v = 50 and x = 1/4 and normalize it so that it goes to zero at a; = 0.) at a given time, and the 
straight lines A, B, C represent the left hand side of 12911 with — "l>|j = 0, k — 1 — e'" , k ~ 1 + e~", respectively, for k > I. 
The arrow indicates how the lines move as increases. The potential has minimum at s = 0. 




FIG. 4: Typical form of the effective potential V at a given time t. The horizontal axis is q. The global minimum is at q 7^ 
and g = is a local maximum. 



and X is still determined by Eq. 129|l , whose solution is obtained by graphical method shown in Fig. 3 where the form 
of the potential is different from that in Fig. 1. With these we schematically plot the potential in q in Fig. 4, which 
shows the presence of the symmetry breaking clearly. A few comments are in order. First, the potential is plotted 
as a function of q. Since the minimum value of q corresponds to g = 0, the form of the potential can roughly be 
interpreted as a potential for q if > 0. On the other hand, if < 0, the relation between q and q changes in time. 
Therefore, we should interpret the potential only at a given time. Second, the restriction k > in the potential H38|) 
disappears in the present case. Even for negative k, the solution at x = 1 corresponds to the smaller potential value 
of the two positive solution of (|29|l for large enough v. 

All of the discussions in this section are the zero temperature results. The inclusion of non-zero initial temperature 
does not leads to a critical complication in the analysis. The only difference is that all the integrals in the calculations 

should contains the coth — — factor. This results in the change of I — > Tcoth — — for the unstable modes and 

2 ^2 
the V in Eq. H38fl becomes temperature dependent. However, these cannot alter the global behavior of the effective 

potential. 
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VII. EQUAL TIME CORRELATION FUNCTION 



Until now we have calculated the effective action with the precarious renormalization condition. However, the 
correlation function, 



kdk sin kr G'k(i) 



(113) 



1 

1 

27r2r 



/ kdk sin kr G^it) + / kdk sinfcr G^it) + / fcdfc 

-/ m J mf 



l(t) 



where r = |x — y|, is renormalization scheme independent in the sense that the renormalization is related only to 
the stable modes and the dominant contribution to the correlation function comes from the exponential growth of 
the unstable modes. This observation enables one to ignore the last integral in Eq. (|113|l which is related to the two 
point correlation function of stable system, which we ignore in this section compared to the first two integrals which 
contain exponentially increasing factors for large time. 

Let us consider the first integral of Eq. (113) where the dominant part of the integrand is proportional to k'^e^^'^'^^\ 
which is peaked around k where k~^ = dt'uj^^{t'). Then, the correlation function H113() can be approximated by 
the contribution of the unstable mode G'^, 

1 



G^(x,y;i) 







27r2 ^ 

px'3/2 



edk ^ Gl{t) 



kr 



1 

8^ 



kr 



(114) 



8 



n 



— ikr 



kuJk 2ir k*0Jk* 2ir 



where k = k + i — . For large i, k/m becomes very small, and the coefficient of /c^e^^''*^*'' should be smooth for small k. 

Another consequence of the small k is that the Gaussian integral well approximates G^ for most range of < t < T, 
until fh > k. This is why we approximate A;2g2ek(t) g^g g, Gaussian function (the steepest descent method) factoring 

sin kr , , 

out the smooth function — in the integrand. 

kr ^ 

Let us consider the r ^ limit of this function. 



limG^(x,y;i) = — 



3/2 ^2ej(i) 



2tt 



n 



+ .k 



2tt. 



1 ^n+^ 

— lim 3 — 

2 r^O rKLOii 



^G^(x,x;t) = ^G(0;i), 



(115) 



n 



1 d n 



where in the last equality we use lim 5 

r— >0 

G'^(x, x; t) = i|i(.|<m (^)- Therefore in general we write 



r->o mujf^ 27 dk \ kujj. 



k'^&f, 



At r = 0, this should correspond to 



G(x,y;0<t<r) ^ lG^(0;i)^^ 



exp 



_^ikr 



n 



-ikr 



rKLU^{t) 



rK*u}f,*{t) 



(116) 



where most of time-dependencies are given by G'^(0;t) in Eq. ^1^. Note that l/fc^ ~ 2 dt'Cii^^t') mcreases as 
time. The exponent Ai^r^/S determines the correlation length. 
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^D{t) = 2V2 


f dt'c,i\t') 

Jo 


~ 2V2 


f dt'rh-\t') 
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which gives the Cahn- Allen relation for constant m^, time lagged deformed relation 4 



for linear mass-squared, ni?{t) 



i (m/ - ^7n] ~ at) 



(117) 
-,1/2 



at, which was discovered in Ref. [T5l |. and another deformed relation 



20 



2 

— tan 

a 



at 



aH 



2+2 



1/2 



for quadratic mass-squared m^it) = 



a^t^. Note that the last two deformed 



and \ —, respectively. This is because the in- 
a V a 



relations have the maxima of the correlation length given by 

stability ends to increase at the time where m^(t) = 0. An interesting feature of the last relation is that the maximum 
correlation length does not depend on the initial mass-squared m^, rather, it depends only on the acceleration of the 
change of the mass-squared. Since we do not know the exact time-dependence of (t) due to the limitation of WKB 
method mentioned in Sec. V, we cannot predict the exact form of the correlation length. 

Next, let us consider the second integral of Eq. 1113|) . The dominant part of the integral is proportional to Pe^^''^*'=\ 

— L J . The approximation of this integral is illustrated in Appendix B. Since the integral 
_ at J 

is maximized at fc ^ to/, we can ignore this term until fhit) ^ k ^ 0. Therefore, during the process of instability 
growing, < t < T, we can ignore the second integral. In this sense, this integral becomes important only after t > T 
since we have no term in this region. The dominant contribution to (^(x, y; t) comes from the second integral for 
t > T and it gives 



kdk sinfcr G^{t) 



1 

2^ 



^/ - o 



sin kr 



\2TTatJ 
\2TratJ 
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exp 



2uj]i kr 
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[l + sin2(?k(t)] 
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, r \ e"" I - .k^r 

q \ k -\- 1 q \ k — i 
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,2efc(t) 



-2© U[fc + z^e(t,-r)] 



J[2eT (t)-kr] \ ^ 



^i[29r(t)+fcr]' 



2r 



(118) 



-ikr 



2ir 



where g{k) 



2kuJu 



and $(t,r) is given in Eq. (jB10|) . Because of the exponential decaying factor, most terms of 



G'(x, y;t) decrease to zero for large r and t. However, a non-trivial exception exists which signals the end of the 
spinodal line along the line 8(t, —r) — since, on this line the exponential decaying factor e^^® of the last line 

in Eq. (|118|l disappears, which gives the non-trivial long time behavior of the correlation function after symmetry 
breaking. 



G(x,y;t) 



\2TTatJ 



3/2 



exp 



k^r\ e''^'" /- k^r 

i q { k — i 

A J 2ir ^ \ 4 



.^.3/2 fVl^2es(t),-2e^(t,-r) 

\2natJ 
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cos [29^ (t) - kr] 



kr 



ikr 



2ir 



(119) 



The first term is the usually expected time-independent correlation, which corresponds to a domain formation and 
growth. The second term of Eq. (|119|l is a new time-dependent correlation, which propagates even after the end of 

the phase transition, and travels along the line kr = — J- f — — — . This term may be interpreted as a propagating 



domain. 



VIII. SUMMARY AND DISCUSSIONS 



In this paper we have elaborated the extension of the double Gaussian type wave-functional approximation to study 
the non-equilibrium quantum dynamics of self-interacting scalar field system in (3-f 1) dimensions. The time-dependent 
second order phase transition is a classic example of this type of problems. These systems are characterized by time- 
dependent coupling parameters and their true nonequilibrium evolution deviates significantly from the equilibrium 
one when their coupling parameters differ greatly from their initial values. In this case the systems evolve completely 
out of equilibrium. To understand the dynamical aspects of these processes there have been developed many different 
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methods such as the closed-time path integral method, sometimes in conjunction with the large TV expansion, mean- 
field, Hartrce-Fock, and the Liouville-von Neumann methods. The present paper uses the functional-Schrodinger 
method and WKB approximation in connection with the Liouville-von Neumann approach to obtain the effective 
action of explicitly time-dependent system undergoing phase transition. 

By applying this method to the scalar theory undergoing second order phase transition, we have found the 
effective equations of motion and effective action of the system in terms of the correlations of the field. This equation 
of motion was used to obtain the effective Hamiltonian of the system written in terms of the mode solutions, which 
is the starting point of the present method. To have better approximation of the theory, we first apply the method 
to the case of constant mass-squared including the tachyonic modes. To do this, we renormalized the effective action 
by employing the standard precarious renormalization scheme and show that the effective potential has no symmetry 
breaking even if one includes the tachyonic modes. The inclusion of the tachyonic modes only results in the presence 
of a new metastable state at the zero mass region. Since what we would like to describe is the system with time- 
dependent mass-squared which evolves in time self-consistently, we added the exact time dependent mode solutions 
under a prescribed change of the mass-squared. By extracting the limiting behaviors of this mode solutions we get the 
WKB approximated solutions of the general time-dependent self-interacting scalar field system with instantaneous 
quenching at t = 0. The full renormalization of the equation of motion and the effective action was performed to 
obtain the finite dynamical evolution of the zero-mode and the mass-squared. By analyzing the time-independent part 
of the effective potential we show that the vacuum structure is changed from that of the system without quenching. 
It is shown that the effective potential indicates the existence of the second order phase transition, which have a new 
vacuum at the non-zero position of the field expectation value and the position of zero expectation value corresponds 
to the unstable equilibrium. 

Due to the limitation of the WKB approximation, we cannot predict the time-dependence of mass-squared if it is 

negative. On the other hand, the effective potential and the spatial correlation function arc calculable. For these 
calculations, we develop a large instability approximation which takes into account the exponential increase of the 
unstable modes dominating the Hamiltonian after renormalization. We computed the spatial correlation function 
with this approximation. Before the phase transition ends, the correlation is dominated by the unstable modes with 

very low frequency determined by the formula k = dt'rh ^(t') , which decreases as time. The inverse of k 

corresponds to the correlation length of the system. We show that the classical Cahn- Allen relation holds only when 
the mass-squared is a negative constant and we have shown that the correlation length depends in general on the 
time-dependence of the mass-squared. An interesting addition in this paper is that the maximum correlation length 
exists depending on the time-dependence of the mass-squared. Especially, if the mass-squared increases quadratically 
in time, the maximal correlation length is independent of the initial mass-squared. Once the phase transition process 
ends, the spatial correlation is dominated by the stable modes with low frequency. Besides the usual static correlation 
corresponding to the formation and growth of domains, it is shown that there exists a travelling correlations. This 

/ TTl f 1 

correlation starts to expand spherically when the phase transition ends with the travelling velocity v ~ 2W — ^ — -— . 

V at 

2V2 UJu - 

The spatial size of this correlation is - . Therefore, this correlation is separated from the static one in time — —at 
_ k ruf 

where t T. The presence of this correlation can be tested on experiment. 
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APPENDIX A: LONG TIME APPROXIMATION 



The function f{y) of Eq. (36) is given by 



f{y) = 2TOt(l-2/)i/2+lny-iln(l-y)+lnm, < y < 1. 



Its first and second derivatives become 

f'iy) = 

f"{y) = 



1 



2y(i - y) 

fht 



2-y- 2mty^Jl - y 



+ ■ 



2(i_y)3/2 y2 2{l-yf 



(Al) 

(A2) 

(A3) 



The function g{y) = 2 — y — 2mty^^/\ — y in the parenthesis of /' has 2 roots in the defined range of y since g(0) = 2. 
g(l) = 1, and g(l/2) < for large fat. Since g{y) decreases at y = 0, smaller solution to this corresponds to the 

maximum value of f{y). In the large t approximation, the smaller root is y = Then, /(y) can be series expanded 
as 
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1 \ 



1 



V 



2 \ mt ) mi ^ 

Then, the integral for slowly varying function g(k) becomes 
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rht 
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{mtf 
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dkg{k)e 



~ g{k)e 



dk exp 



m \ \ t 
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APPENDIX B: LARGE INSTABILITY APPROXIMATION 



The first step of the large instability approximation is to approximate f{y) of Eq. (99) up to quadratic part around 
the maximum point of it: 



f{y) 



= 2m/ ^ ^u{t') - ydt' + Iny + 21nm/ = f{y) -^iy-yf + ' 



where vit) = m?(t)/m'j, {y,t = t^) is the point where f{y) is maximum determined by 



1 f dt 

y Jo ^/iy{t) - y 



and 7 is the negative of the curvature of / at that point: 



j=-dlf{y)\ _ = l + !^^lim 



dt' 



Ht')-yf/^ i,{t)^^4ty^ 



(Bl) 



(B2) 



(B3) 



which is 0{m'jP). Let us approximately estimate the size of y. Let the rate of change of the mass squared m^{t = 



r) = 0. We can safely assume that the acceleration 



d^m^it) . 



dt^ 



is non-negative during t < t < T, which can be easily 



understood from the potential in Fig. 1 , where all unstable modes tend to increase the mass-squared if it is negative. 
Therefore, with ^{t) = fh'^{t)/m'j, ^{t) — y satisfies (1 — y)(l — t/tk) < v{t) — y < 1 — 2/ since its initial value 
i'{t = r) = 1 and ^{tk) — y = 0. With this inequality we get 



is 



tk 



y < I dt^/v(t) - y < tk\/l - y 



< / ^=^= < 



dt 



2tk 



Jo ^v{t) - y Jo ^y{l-y){l-t/tk) 



(B4) 
(B5) 



23 



which leads to, 



with 1 < a <2. With this result, 7 ~ a'^mj-P. 
Since this 7 is large, we can approximate 



2 - 72 



at ' 



TTif \ \ at 



Series expanding J^* uj\i^{t')dt' to first order in (fc — k) we have. 



(B6) 
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u;i,{t')dt' = / Wk(i')'^^'+ / [cJk{t') - UJk{t')]dt' + / uJ-k{t')dt 



—j{k-k) + j_ u;-,{t')dt 



(B8) 



t - k - 
where we ignore /. uji^{t')dt' since it is 0{k) for all time and we use tDk — Wj, H {k — k) + ■ ■ ■ in the second equality. 

Using this equation we get 

— {k-kf±2ie-{t)±ikr ~ — \k-k±ikQ{t,r)f + 2Q'^{t,r)±i\kr + 2ei{t)\ , (B9) 
m/ TO/ L fe J 

where Q{t,r) is dimensionless function. 



e(i,r) 



TO/ /■* di' fcr 



2atJ-, u-k{t') 4 



(BIO) 



Note that 8(t) = 0(t,r = 0) becomes significant only after t — t > t. With these we can write some general form of 
the approximation of the integral, 
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g [k±ikQ{t,r)] e 
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APPENDIX C: ROUGH ESTIMATION OF G INTEGRALS 



The integral Q{m,t) in Eq. (104) 

G^(t)-/o(TO2) 

|k|>m/ 



9^ 



|k|>, 
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k=0 2wk 



|k|>r 



2lji 



■cos26'+(i). 



(CI) 



consists of three terms. The final term is unimportant to our calculation and decreases as time so we ignore it. The 
second integral is exactly integrable to give 



Z"™/ 1 mj 
Jk=o 2wk Stt^ 



1 + ^ 



(C2) 



The first integral of Eq. ljCl|l cannot be exactly integrable. Since the argument of the integral decreases rapidly, 
see ()89f) . we series expand the integrand around k — ruf: 
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Since gj^ — 1 is always positive definite, the above series breakdown at the point, k' = ruf + 
series becomes zero. Therefore, we estimate the first integral by 
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In summary. 
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where g{xi. xj) 



1 



2407r2 \xf 
Similarly, we may obtain 
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where g'{xi,Xf) is a positive number dependent on Xi and only. Since we only need its large m dependence, we 
do not explicitly write g' . 

Other integrals which we need in calculating in Eq. (109) are given by 
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